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Abstract 

We present a graphical representation of the supersymmetry and a 
C-program for the graphical calculation. Calculation is demonstrated 
for 4D Wess-Zumino model and for Super QED. The chiral operators 
are graphically expressed in an illuminating way. The tedious part of 
SUSY calculation, due to manipulating chiral suffixes, reduces consid- 
erably. The application is diverse. 
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1 Introduction 

The supersymmetry is the symmetry between fermions and bosons. It was 
introduced in the mid 70's. At present the experiment does not yet confirm 
the symmetry, but everybody accepts its importance in nature and expects 
fruitful results in the future developement. The requirement of such a high 
symmetry costs a sophisticated structure which makes its dynamical analysis 
difficult. In this circumstance, we propose a calculational technique which 
utilizes the graphical representation of SUSY. The representation was pro- 
posed in [HO]. 1 The spinor is represented as a slanted line with a direction. 

*On leave of absence from Laboratory of Physics, School of Food and Nutritional 
Sciences, University of Shizuoka. Untill 31 March, 2006. 
E-mail address: ichinose@u-shizuoka-ken.ac.jp 

J An improved version of Ref.pp has recently appeared as Ref.jS]. 
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Its chirality is represented by the way the line is drawn. The introductory 
explantion is given in the text. The advantage of the graph expression is the 
use of the graph indices. Every independent graph, which corresponds to a 
unique term in the ordinary calculation, is classified by a set of graph indices. 
Hence the main efforts of programinng is devoted to find good graph indices 
and to count them. SUSY calculation generally is not a simple algebraic or 
combinatoric or analytical one. It involves the vast branch of mathematics 
including Grassmann algebra. The delicate property of chirality is produced 
in this environment. Hence it seems that the ordinary(popular) programs, 
such as Mathematica, Maple, REDUCE, MAXIMA, form, do not work. It 
requires a more fundamental language. We take C-language and present a 
first-step program. Future development to the level of the previously cited 
ordinary programs is much expected. 

The notation in the text is based on the standard textbook by Wess and 
Bagger |3] . 

2 Spinors in SUSY: Graphical Representation 
and Storage Form 

Weyl spinors have the SU(2)lxSU(2)^ structure. The chiral suffix a, appear- 
ing in tp a or ip a , represents (fundamental representation, doublet representa- 
tion) SU(2)l and the anti-chiral suffix a, appearing in ip a or ipa, represents 
SU(2)#. The raising and lowering of suffixes are done by the antisymmetric 
tensors e alS and e a p. 

«*) = ( _°i I ) • <«w> = ( ? "o 1 ) • = s " • 

= , <fc= t4( # . (1) 

They are graphically expressed by Fig.l. We encode them as follows. We use 
2 dimensional array with the size 2x2. The four chiral spinors are stored in 
C-program as the array psi[ ][ ]. 

Fermionic Fields 

(1) (Weyl) Spinor [Symbol: p ; Dimension: M 3//2 ] 



1p a 
















psi[0,0] = 


=a 


psi[0,0] = 


=empty 


psi[0,0] = 


=empty 


psi[0,0] = 


=empty 


psi[0,l]= 


=empty 


psi[0,l] = 


=a 


psi[0,l] = 


=empty 


psi[0,l] = 


=empty 


psi[l,0]= 


=empty 


psi[l,0] = 


=empty 


psi[l,0]= 


=empty 


psi[l,0] = 


=a 


psi[l,l]= 


=empty 


psi[l,l] = 


=empty 


psi[l,l]= 


=d 


psi[l,l] = 


=empty 
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Figure 1: Weyl fermions. 



The first column takes two numbers and 1; expresses a 'crural' operator 
if), while 1 expresses an 'anti-chiral' operator if). The second column also takes 
the two numbers; expresses an 'up' suffix, while 1 expresses an 'down' one. 

Note: Chiral spinor suffixes a, (3, - ■ ■ are expressed , in the present de- 
program, by positive odd number integers 1,3,- ■ ■ , while anti-chiral ones 
d,P,--- are by positive even number integers 2,4,- • ■ . This convension (dis- 
criminative use of even and odd integers) is, at this stage, rather redundant 
in the sense that the chirality (if) or if)) can be read by the first element 
number of the array psi[2][2] for a non-empty data, i.e. for if) (psi[0] [*]=«) 
and 1 for if) (psi[l][*]=d). (The situation is the same for some other spinors 
d m if>, 9, ■ ■ ■ . See the later discription. ) The convension, however, will soon 
become important to discriminate the chirality of the spinor matrices; a and 
a. 

Note: 'empty' is expressed by a default number (, for example, 99) in the 
program. 

(2) First derivative of spinor [Symbol: q ; Dimension: M 5 / 2 ] 
The first derivative of the spinor is graphically expressed by the upper graph 
of Fig. 2. It is stored by one 2x2 array dps[ ] ] and one variable dpsv. 
d m if) a d m if) a d m ipa d m ip a 

dps[0,0]=a dps[0,0]=empty dps[0,0]=empty dps[0,0]=empty 

dps[0,l]=empty dps[0,l]=a dps[0,l]=empty dps[0,l]=empty 

dps[l,0]=empty dps[l,0]=empty dps[l,0]=empty dps[l,0]=o; 
dps[l,l]=empty dps[l,l]=empty dps[l,l]=d dps[l,l]=empty 
dpsv=m dpsv=m dpsv=m dpsv=m 

Here the vector suffix expresses the Lorents suffix of a differetial operator. 

Note: The vector suffixes m, n, ■ ■ • are expressed by 51, 52, ■ ■ • in the 
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Figure 2: Derivatives of fermions. 



present program. 



(3) Sigma Matrix [Symbol: s ; Dimension: M° ] 
Sigma matrices a m , a m are graphically expressed in Fig. 3. 
They are stored as the 2x2 array si[ ][ ]. 



a m ■ 

" aa 



si[0,0]=empty si[0,0]=Q! 
si[0,l]=a si[0,l]=empty 
si[l,0]=empty si[l,0]=a; 
si[l,l]=d si[l,l]=empty 
siv=m siv=m 

Note: The use of even (a) and odd (a) integers makes an important role 
here. The arrangement of spinor suffixes, that is (left-side suffix, right-side 
suffix) = (odd, even) or (even, odd), makes us clear the difference between a 
and a. We should, however, have the relation 



—m m 

aa ad 



in mind. Hence the above 2quantities are equivelently expressed as 



" aa 



si[0,0]=empty si[0,0]=o; 
si[0,l]=d si[0,l]=empty 
si[f ,0]=empty si[f,0]=(i 
si[l,l]=a si[f ,1 ]=empty 

siv=m siv=m 

This ambiguity does not cause any problem because we keep a rule: 
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Figure 3: Elements of SL(2,C) a-matrices. (c m ) Q/ 3 and (<j m ) alS are the 
standard form. 



In this program, we use only a (not use a r 



• a rn is used only for the graphical explanation. 

(4) Superspace coordinate [Symbol: t ; Dimension: M" 1 / 2 ] 

The superspace coordinate 9 a is exprssed in the same way as the spinor ip a . 







9 a 




9a 








th[0,0]= 


=a 


th[0,0]= 


=empty 


th[0,0]= 


=empty 


th[0,0]= 


=empty 


th[0,l]= 


=empty 


th[0,l]= 


=a 


th[0,l]= 


=empty 


th[0,l]= 


=empty 


th[l,0] = 


=empty 


th[l,0]= 


=empty 


th[l,0]= 


=empty 


th[l,0]= 


=a 


th[l,l] = 


=empty 


th[l,l]= 


=empty 


th[l,l]= 


=a 


th[l,l]= 


=empty 



They are graphically expressed by Fig. 4. 
(5) Gagino [Symbol: 1 ; Dimension: M 3 / 2 ] 

The photino A" is exprssed in the same way as the spinor tfj a . We take the 
2x2 array laj ][ ]. 
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Figure 4: The graphical representation for the spinor coordinates in the 
superspace: 9 a ,9 a ,6a and 6 a . 



la[0,0]=a 
la[0,l]=empty 
la[l,0]=empty 
la[l,l]=empty 



la[0,0]=empty 
la[0,l]=a 
la[l,0]=empty 
la[l,l]=empty 



Aq, 

la[0,0]=empty 
la[0,l]=empty 
la[l,0]=empty 
la[l,l]=d 



la[0,0]=empty 
la[0,l]=empty 
la[l,0]=d 
la[l,l]=empty 



(6) The first derivative of gaugino [Symbol: m ; Dimension: M 5 / 2 ] 

The first derivative of the photino is expressed as the 2x2 array dl[ ] ] and 

the variable dlv. 

d m X a d m )\ a d m Xa d m X a 

dl[0,0]=a dl[0,0]=empty dl[0,0]=empty dl[0,0]=empty 

dl[0,l]=empty dl[0,l]=a dl[0,l]=empty dl[0,l]=empty 

dl[l,0]=empty dl[l,0]=empty dl[l,0]=empty dl[l,0]=d 
dl[l,l]=empty dl[l,l]=empty dl[l,l]=d dl[l,l]=empty 
dlv=m dlv=m dlv=m dlv=m 



Bosonic Fields 

(7) Complex scalar [Symbol: A ; Dimension: M 1 ] 

The complex scalar field A is expressed by one dimensional array A[ ] with 

2 elements. 

A A* 

A[0]=l(exist) A[0]=empty 
A[l]=empty A[l]=l (exist) 

where the element-numbers 0,1 correspond to ^4(chiral) or A*(anti-chiral), 
respectively. 
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(8) The first derivative of the complex scalar [Symbol: B ; Dimension: M 2 ] 
The first derivative of A and A* are expressed as the one dimensional array 
B[ ] with 2 elements. 

B[0]=m B[0]=empty 
B[l]=empty B[l]=m 

(9) Vector field [Symbol: v ; Dimension: M 1 ] 

The vector field (photon) v m is expressed by one variable v. 

v m 
v=m 

v=empty (non-exist) 

The lowest expression is taken when the vector field does not appear. 

(10) The first derivative of the vector field [Symbol: w ; Dimension: M 2 
] 

The first derivative of v m is expressed by two variables dv and dvv. 

dv=m 
dvv=n 



(11) The Dalemberian derivative of A and A* [Symbol: C ; Dimension: M 3 ] 
The Dalemberian derivative of A and A* are expressed as the one dimen- 
sional array C[ ] with 2 elements. 
d m d m A d m d m A* 
C[0]=1 (exist) C[0]=empty 
C[l]=empty C[l]=l(exist) 



(12) Auxiliary fields [Symbol: F ; Dimension: M 2 ] 

The auxiliary fields F and F* are expressed as 

F F* 

F[0]=1 (exist) F[0]=empty 

F[l]=empty F[l]=l (exist) 



(13) Real auxilary fields [Symbol: D ; Dimension: M 2 ] 

The auxiliary field D (real scalar), which appears in the vector multiplet, is 

expressed as 
D 

D=l (exist) 
D=empty (non-exist) 



Term and Component 

Terms in the SUSY calculation are stored as an ordered set of the above 
quantities. The following examples appear in the intermidiate stage of eval- 
uating W a W a where W a is the field strength superfield. (See App.B.) 

Example 1 

The term a l ^a m ^a s5 -ay 8^8^= -a^a^a^a^d^d^ is stored, 
in the computer, as follows. 

type[c=0]=s type[c=l]=s type[c=2]=s type[c=3]=s 
si[c=0,0,l]=l si[c=l,0,0]=3 si[c=2,0,0]=l si[c=3,0,l]=3 
si[c=0,l,l]=2 si[c=l,l,0]=2 si[c=2,l,l]=4 si[c=3,l,0]=4 ^ed*" 
siv[c=0]=51 siv[c=l]=52 siv[c=2]=53 siv[c=3]=54 below 



type[c=4]=w type[c=5]=w 
dv[c=4]=52 dv[c=5]=54 
dvv[c=4]=51 dvv[c=5]=53 

(8 -> 1, a -> 2, -> 3, 7 -> 4; 51, m -> 52, s -> 53, n -> 54) 
Besides the above ones, the number of components (compno=6) and an over- 
all weight (complex) are necessary to characterize a term. One additional 
coloumn, specified by the variable c, appears. The variable c controls the 
order of every component. This ordering is important in the calculation in- 
volving Grassmannian quantities. In the above storing form, each component 
starts with specifying the type: s, w, • • • . The one dimensional array type[ ] 
is used for the purpose. This term graphically appears in the output (before 
further reduction) as shown in Fig. 5. 




Figure 5: The graphical expression of — cr l Sdl a malS a^a^ 6 div m d s v n . 

In Fig. 5, we see all dummy suffixes disappear and the advantage of the 
graphical expression is manifest. The chirality can be read from the shape 
of the directed-line graph. 

Example 2 

The term o\a m 8 a Ddiv m = —a l S6c a ma5 Ddiv m is stored, in the computer, as 
follows. 
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type[c=0]=s type[c=l]=s 
si[c=0,0,0]=l si[c=l,0,l]=l 
si[c=0,l,l]=4 si[c=l,l,0]=4 type[c=2]=D 
siv[c=0]=53 siv[c=l]=54 D[c=2]=l 
(<J->l,d->4; /^53,m^54) 

In this case, the number of components is 4 (compno=4) and there is a 
weight. This one graphically appears in the output (before further reduction) 
as shown in Fig. 6. 



type[c=3]=w 

dv[c=3]=54 

dvv[c=3]=53 




Figure 6: The graphical expression of — a l 56i a ma& Ddiv 



T<Jm- 



These two examples clearly show the advantage of the graph representa- 
tion over the conventional one using indicies. The new representation dis- 
criminates each term not by the suffixes but by the shape of the graph. It 
means the importance of indices of the graph such as the number of the 
chiral-loop. (In the above examples, the number is 1.) 



3 Field name, id-number and Grassmannian prop- 
erty 

One advantage of C-language is the efficient manipulation of characters. The 
control of characters is important for specifying various operators (fields) ap- 
pearing in SUSY theories. As described previously, we assign every operator 
a symbol which is expressed by a character. The appearance of many fields, 
which is a big cause complicating SUSY theories, can be neatly handled by 
the use of the character name. We also introduce the id-number for every 
operator. (See the list below) 

For every operator, we assign the Grassmann number: +1 (commuting) 
or -l(anti-commuting). This assignment is exploited in the process of mov- 
ing a component , within a term, respecting the Grassmannian property of 
operators. 

We list the assignment in TABLE 1 with the dimension of operators. 
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Fid[0]= 


't'; 


DIM[0]=0.0; 




Gra ssm a n n fOl = 

\J J- ( A k J k ) X X X < V X X X X 1 \J 1 


-1- 


Fid[l]= 


V; 


DIM[1]=0.0; 




Gra ssm a nn M 1 — 

V \ J. (XjkJLJ XXXIXXXXX 1 X. 1 


1 ■ 


/* <Chiral */ 










Fid[100]= 


='p'; 


DIM[100] = (float)3/2; 


\JL dOOlIldllll 1UU - 


- i 
j. 


Fid[101]= 


='q'l 


DIM[101] = (float)5/2; 


C-rrassmaTiTi [1 (11 1 = 

V_. ' J. 1 1 k H JXXXCVXXXX XU X 


=-1 


Fid[102]= 


='A'; 


DIM[102]=1.0 




Gra,ssmann[102l : 

\j x cxjkjkjxxx cxjxx hi x \j i 


=1: 


Fid[103]= 


='B'; 


DIM[103]=2.0 




Grassmann [1 031= 


=1: 


Fid[104]= 


='C; 


DIM[104]=3.0 




Grassmann [1 041= 

V ' X CXkJUlllCXll 11 1 X \J J- 1 


=1: 


Fid[105] = 


='F'; 


DIM[105]=2.0 




Grassmann[105] : 


=1: 


'* <Chiral2 */ 










Fid[120]= 


='P'; 


DIM[120] = (float)3/2; 


Grassmann[120]= 


=-1 


Fid[121]= 


='Q'; 


DIM[121]=(float)5/2; 


Grassmann[121]= 


=-1 


Fid[122]= 


='a'; 


DIM[122]=1.0 




Grassmann[122]= 


=1: 


Fid[123]= 


='b'; 


DIM[123]=2.0 




Grassmann[123]= 


=1: 


Fid[124]= 


=V; 


DIM[124]=3.0 




Grassmann[124]= 


=1; 


Fid[125]= 


=T; 


DIM[125]=2.0 




Grassmann[125]= 


=1; 


/* < Vector */ 










Fid[140]: 


=T; 


DIM[140] = (float)3/2; 


Grassmann[140]= 


=-1 


Fid[141]= 


='m'; 


DIM[141]= (float) 5/2; 


Grassmann[141]= 


=-1 


Fid[142]= 


='D'; 


DIM[142]=2.0 




Grassmann[142]= 


=1: 


Fid[143]= 


=V; 


DIM[143]=1.0 




Grassmann[143] : 


=1: 


Fid[144]= 


=V; 


DIM[144]=2.0 




Grassmann[144]= 


=1: 



TABLE 1 Definition of FieldID, Dimension, Grassmann No 
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4 Graph Indices: vpairno, NcpairO, NcpairE, 
closed-chiral-loop-No, GrNum 

SigmaContraction (B),(C),(D) and (E) 

In the process of SUSY calculation, there appear graphs connected by 
directed lines (chiral suffixes contraction) and by (non-directed) dotted lines 
(vector suffixes contraction). We can classify them by some graph indices. 

vpairno The number of vector-suffix contractions. 

NcpairO The number of chiral-suffix contractions. This is equal to the 
number of left-directed wedges. 

NcpairE The number of anti-chiral-suffix contractions. This is equal to the 
number of the right-directed wedges. 

closed-chiral-loop-No The closed-chiral-loop is the case that the directed 
lines , connected by a or a, make a loop. In this case NcpairO=NcpairE. 
The number of closed chiral loops is defined to this index. 

GrNum A group is defined to be a set of cr's or tr's which are connected by 
directed lines. The number of groups is defined to be GrNum. 

In TABLE 2-4, we list the classification of the product of cr's using the 
graph indices defined above. 

These tables clearly show the a-matrices play an important role to connect 
the chiral world and the space-time (Lorentz) world. 



5 Superspace coordinates 

Supersymmetry is most manifestly expressed in the superspace (x m ,9,9). 

9 a = e afi 90 are spinorial coordinates. They satisfy the relations 
graphically shown in Fig. 7. These relations are exploited in the program in 
order to sort the SUSY quantities with respect to the power of 99 and 99. ( 
For further detail, see the subsection thBthBthth of App.A. ) 



6 Treatment of Metrics: e a p, ej* = e^ a = 6%, f] mn 
and the totally anti-symmetric tensor e lmns 

For the totally anti-symmetric tensor e lmns , we introduce one dimensional 
array ep[ ] with 4 components. 
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vpairno 


NcpairO 


NcpairE 


figure 




n 
u 


n 
u 


a a p p 


n 
u 


n 
u 


i 
i 


m - 11 




1 





-m 11 

/v'\ 




1 


1 












a a p p ^ 


1 





1 


a /°VN ="4^ 

a p 




1 









1 


1 





TABLE 2 Classification of the product of 2 sigma matrices (nsi=2). 
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vpairno 


NcpairO 


NcpairE 


figure 










1 m n 
/ CT \ / CT \ / CT \. 

r i a d P p 








1 


1 -m n 

/ CT xy CT \ / G \. 

a By Y 




1 





-1 ra n 

/ a x/ a \ / a \, 

a P Y Y 




1 


1 


closed-chiral-loon No =1 

</ \yj> / \ 








closed-chiral-loop No =U 








1 -m n 

/"V °v \ 

a a 










1 

Y Y 


1 





1 


-4^ /°\, 

Y Y 




1 





n 

-m / a \, 

P y y 




1 


1 


1 

-8 / a \ 

Y Y 



TABLE 3 Classification of the product of 3 sigma matrices (nsi=3). 
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vpairno 


NcpairO 


NcpairE 


figure 











/ \ / \ / \./ \. 

y y a ap p8 S 





1 


A/ a x / a \ /\ 

a P y 7 8 5 


1 





/ 5 V°\. / c \, / G \, 

a P Y y 5 5 


1 


1 


GrNum=2, Division=(2,2) 

/V°\, /\/\ 

a M p 


GrNnm=2 Division=(3 1 ) 

s 

— 1 m — n rr 

/°v"\/°\ / \ 

a a v y 


GrNum=3, Division=(2,l,l) 

1 _ m _ n s 

</^> / V \ 


2 





.. CT x / CT v / CT \ / CT \v 

a M 5 





2 


/ v x/ v \ 

a p 


1 


2 


1 — m n — : ' ; 

GrNum=l, / \/ V/ \/ X 


GrNum 2, / \A / V ^ 


2 


1 


— m n — b l 

GrNum=l, / V V X/ \. 

a p 


1 _m _H s 

GrNum=2, / \/ \/ V/ \ 


2 


2 


GrXum=l, \/\/J> 


1 - m c n 

GrNum 2, < ^/^><^ X 



TABLE 4 Classification of the product of 4 sigma matrices with 
no vector-suffix contraction (nsi=4, vpairno=0). 
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e \\ 9 - -h"" e^e 



a p 



a p 

Figure 7: The graphical rules for the spinor coordinates: 9 a 9 13 
—^e a ^99, 9 a 6p = \t a p06, Q a W = ^e a @99, 9^9^ = 9a m 99a n 9 
-\9999r] mn . 

ep[0]=l 
ep[l]=m 

G ^P] 11 Symbol: e 

ep[3J=s 



This term appears in SecEJand produces topologically important terms such 

as Vl m V lm = e lmnS Vl m V ns . 

As for the metric of the chiral suffix, we do not introduce specific arrays. 
They play a role of raising or lowering suffixes, which can be encoded in 
the upper (0) and lower (1) code in arrays. For the Lorentz metric r] mn , 
we do not need to much care for the discrimination between the upper and 
lower suffixes because of the even-symmetry with respect to the change of 
the Lorentz suffixes (?7 mri = r) nm ). 



7 Sigma matrices 

Let us express important relations valid between products of sigma matrices 
graphically. In Fig. 8, the symmetric combination of a m a n are shown as the 
basic spinor algebra. The antisymmetric combination gives the generators of 
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_m n 



/ \/ \ + m <— > n = -2T| mn 5 a 

p P 



m - n 



/ a \/ G \ + m <~> n =-2^5' 



.P 

'a 

a P 



Figure 8: A graphical formula of the basic spinor algebra. (er m ) a/3 \a n ) g* 



m^n = -2r) mn S^, and {a m )a&(a n ) &f) + m <-> n = -2ry mn 5f. 



the Lorentz group, a nm ,a r 



n — m 



1 1 /v\, 



I —n m ^ 

/V \ -m^nl , (2) 

The "reduction" formulae (from the cubic ex's to the linear one) are ex- 
pressed as in Fig. 9. From Fig. 9, we notice any chain of a's can always be ex- 
pressed by less than three a's. The appearance of the 4th rank anti-symmetric 
tensor e lmns is quite illuminating. The completeness relations are expressed 
as in Fig. 10. The contraction, expressed by arrowed curve in Fig. 10. is the 
matrix trace. 

Fierz identity is graphically shown as follows. 



\ / \ 

U~ /V7 e af3+ \V\ e aP -m~n 



a a p 

-n ^ m =. n 



U /"V 5 / -l~n\\ \ 5 V a '\. -l»m 



Jim 



These relations are important, in the program, to reduce the product of 
sigma matrices. Especially at the place (E) of SigmaContraction, they are 
exploited. As an example, the closed chiral-loop graph in Fig. 5 reduces as 
follows. 

j a. / CT v 

X/ \/ = 2{rf m rf is — rj ln rj ms + rj ls r] mn — ie lmns ) (4) 
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/W ff \ - 

a a 

- / °\ iT + / °\ rf 1 - / <\r\ lm + i e ^ / ^ 

a a a \ a a a 

— 1 m — n 

/ CT v a \/ a \ - 

a a 

/V - / G V m - 

a a a a a a 

Figure 9: Two relations: 1) a l a m a n = —<j l r] rnn + a m r] nl — a n rf m + ie lmns a s , 
2) a'a m a n = -a l r] mn + a m r/ n ' - a n 7] lrn - ze' mns a s . 



m - 11 

/ xy \ = -2-n 

> 

a a ^ p 

Figure 10: Completeness relations: 1) S^{a m ) a( ^{a n )^ = -2r/ mn , 2) 

(a m ) m (ff m )^ = -2<y£<Sj. The contraction using 52 is the matrix trace. The 
relation 1) is also obtained from Fig. 8. 
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Its complex conjugate one is given by 

= 2{rf m rf ls - rf n T] ms + rf s r] mn + ie lmns ) (5) 

The above 2 terms appear in the calculation of W a W a and WaW a respec- 
tively. (W a and Wa are superfields of field strength. ) 

8 Superfield and Structure of Input Data 

The transformation between the superfield expression and the fields-components 
expression is an important subject of SUSY theories. For the purpose, we 
do the calculation of in the App.B. In this case, the input data is 

taken from the content of the superfield. The first superfield <3>t(we assign 
the superfield-number sf=0) and the second one $(sf=l) have 6 terms (we 
assign the term-number t=0,l,- • • ,5) within each superfield. The input data, 
written in App.B, can be read from the following expression. 



m 

& = -i \ / a \/ d m A*+ K 



x / -d 2 A* + 2 \ / 



v - 4 



+i \ / + ± / F* + A* (6) 

This data of is stored as 

weight [sf=0,t=0]=0+i(-l) weight[sf=0,t=l]=l+i(0) 

type[sf=0,t=0,c=0]= t type[sf=0,t=l,c=0]= t 

th[sf=0,t=0,c=0,0,0]=l th[sf=0,t=l,c=0,0,0]=l 

type[sf=0,t=0,c=l]= s type[sf=0,t=l,c=l]= t 

si[sf=0,t=0,c=l,0,l]=l th[sf=0,t=l,c=l,0,l]=l 

si[sf=0,t=0,c=l,l,l]=2 type[sf=0,t=l,c=2]= t 

siv[sf=0,t=0,c=l]=51 th[sf=0,t=l,c=2,l,l]=2 

type[sf=0,t=0,c=2]= t type[sf=0,t=l,c=3]= t 

th[sf=0,t=0,c=2,l,0]=2 th[sf=0,t=l,c=3,l,0]=2 

type[sf=0,t=0,c=3]= B type[sf=0,t=l,c=4]= C 

B[sf=0,t=0,c=3,l]=51 C[sf=0,t=l,c=4,l]=l 

Here we notice two additional coloumns, specified by sf and t, appear. They 
specify the order of the supersields and the order of the terms within each 
superfield respectively. 
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The superfield $ is graphically shown as 

$ = i \y G \/ d m A+ \/ - A d 2 A + 2 s/ 



This data of $ is stored as 
weight [sf=l,t=0]=0+i(l) weight[sf=l,t=l]=l+i(0) 

type[sf=l,t=0,c=0]= t type[sf=l,t=l,c=0]= t 

th[sf=l,t=0,c=0,0,0]=5 th[sf=l,t=l,c=0,0,0]=5 

type[sf=l,t=0,c=l]= s type[sf=l,t=l,c=l]= t 

si[sf=l,t=0,c=l,0,l]=5 th[sf=l,t=l,c=l,0,l]=5 

si[sf=l,t=0,c=l,l,l]=6 type[sf=l,t=l,c=2]= t 

siv[sf=l,t=0,c=l]=52 th[sf=l,t=l,c=2,l,l]=6 

type[sf=l,t=0,c=2]= t type[sf=l,t=l,c=3]= t 

th[sf=l,t=0,c=2,l,0]=6 th[sf=l,t=l,c=3,l,0]=6 

type[sf=l,t=0,c=3]= B type[sf=l,t=l,c=4]= C 

B[sf=l,t=0,c=3,0]=52 C[sf=l,t=l,c=4,0]=l 

The calculation of <3>f<3> leads to the Wess-Zumino Lagrangian. (See 
App.B) 

We also do the calculation of W a W a in App.B where W a is the field 
strength superfield. They are expressed as follows. 

W a = -i A + Q / D-i /VV \n mn+ V / Gm XX 
a a or v 2 v a 

This data of W a is stored as 

weight[sf=0,t=l]=l+i(0) 

type[sf=0,t=l,c=0]= t 
weight[sf=0,t=0]=0+i(-l) th[sf=0,t=l,c=0,0,l]=l 
type[sf=0,t=0,c=0]= 1 type[sf=0,t=l,c=l]= D 
la[sf=0,t=0,c=0,0,l]=l D[sf=0,t=l,c=l]=l 

W a is expressed as 

w = -i V + e\ D-i Xfxy 9 "\/ \* mn+ \J \fxy„ 

This data of W a is stored as 
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weight[sf=l,t=l]=l+i(0) 
type[sf=l,t=l,c=0]= t 
weight [sf=l,t=0]=0+i(-l) th[sf=l,t=l,c=0,0,0]=l 
type[sf=l,t=0,c=0]= 1 type[sf=l,t=l,c=l]= D 
la[sf=l,t=0,c=0,0,0]=l D[sf=l,t=l,c=l]=l 
The calculation of W a W a leads to the Super Electromagnetism Lagrangian.(See 
App.B) 

9 Conclusion 

In the history of the quantum field theory, new techniques have produced 
physically important results. The regularization techniques are such exam- 
ples. The dimensional regularization by 'tHooft and Veltman[H] produced 
important results on the renormalization group property of Yang-Mills the- 
ory and many scattering amplitude calculations. The lattice regularization 
in the gauge theory revealed non-perturbative features of hadron physics. In 
this case, the computor technique of numerical calculation is essential. As for 
the computer algebraic one, we recall the calculation of 2-loop on-shell coun- 
terterms of pure Einstein gravity [12 El- A new technique is equally important 
as a new idea. 

The SUSY theory is beautifully constructed respecting the symmetry be- 
tween bosons and fermions, but the attractiveness is practically much reduced 
by its complicated structure: many fields, chiral properties, Grassmannian 
algebra, etc. The present approach intends to improve the situation by a 
computer program which makes use of the graphical technique. (This ap- 
proach is taken in Ref. |H| for the calculation of product of SO(N) tensors. It 
was applied to various anomaly calculations. ) 

The present program should be much more improved. Here we cite the 
prospective final goal. 

1. It can do the transformation between the superfield expression and the 
component expression. 

2. It can do the SUSY trnasformation of various quantities. In particular 
it can confirm the SUSY-invariance of the Lagrangian in the graphical 
way and give the final total divergence. 

3. It can do algebraic SUSY calculation involving D a , D a , Q a and Q a . 

The item 1 above has been demostrated in the present paper. 

It is impossible to deal with all SUSY calculations. This is simply be- 
cause which fields appear and which dimensional quantities are calculated 
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depend on each problem. If we obatin a list of (graph) indices which classify 
all phsical quantities (operators) appearing in the output, then the present 
program works (by adding new lines for the new problem). To deal with 
such a case, we add the appendix B where the program flow is explained. It 
intends to help the reader to read the original source code. 

10 Appendix A: Skech of Programming Flow 

2 

2 For simplicity, we omit the vector multiplet components. 
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10.1 Input of Data 

[Initialization of Basic ComponentsJ 

sf=0,l,- • • maxsfieldno-1; t=0,l,- • • maxtermno-1; c=0,l,- • • maxcompno-1; i,j,k=0,l 

type[sf][t][c]='x' 

th[sf][t][c][i][j]=99 

psi[sf][t][c][i][j]=99 

dps[sf][t][c][i][j]=99 

dpsv[sf][t][c]=99 

si[sf][t][c][i][j]=99 

siv[sf][t][c]=99 

A[sf][t][c][k]=99 

B[sf][t][c][k]=99 

C[sf][t][c][k]=99 

F[sf][t][c][k]=99 

siieldno=2 

for (sf=0;sf<sfieldno;sf++) 
{ /* <sf-running */ 
printf(" SUPER FIELD NO" ,sf) 
fscanf (& termno) — > printf, notermfsf] 

for (t=0;t<termno;t++) 
{ /* <t-running */ 
printf ( " term= " ,t ) 



[Initialization of Graph Indices] 

dif=0, nth=thn=thbarn=0, npsi=psin=psibarn=0, nsi=sin=sibarn=0 
nA=An=Abarn=0, nB=Bn=Bbarn=0, nC=Cn=Cbarn=0, nF=Fn=Fbarn=0 



fscanf (&weight[sf][t][2]) — > printf 
fscanf (&compno) — > printf, nocomp[sf][t] 

for (c=0;c<lcompno;c++) 
{ /* <c-running */ 

printf(c) 

fscanf (&type[sf 



[t][c]) — > printf 



SWITCH 
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[Install of SuperField Data] 
case Y /* Spinor Coordinate 9, 9 */ 
nth— > nth+1 printf(nth) 
fscanf(&lr,&ud, fespsuf) 

lr=0 thn — > thn+1 printf(thn) 
lr=l thbarn — > thbarn+1 printf(thbarn) 



switch(lr) 



printf (spsuf) 
th[sf][t][c][lr][ud]=spsuf 



case 'p' /* Spinor if),v/j */ 

npsi— > npsi+1 printf(npsi) 
fscanf(&lr,&ud, &spsuf) 



switch(lr)-> { lr= ° pSin ~* psin+1 P rintf (P sin ) 

^ ' \ lr=l psibarn — > psibarn+1 printf(psibarn) 



printf (spsuf) 

psi [sf ] [t] [c] [lr] [ud] =spsuf 



case 'q' /* der Spinor dip, d$ */ 

npsi— > npsi+1, dif— >dif+l printf (npsi) 
fscanf(&lr,&ud, &spsuf) 

lr=0 psin — > psin+1 printf(psin) 
lr=l psibarn — > psibarn+1 printf(psibarn) 



switch(lr) 



printf (spsuf) 

dps [sf ] [t] [c] [lr] [ud] =spsuf 

fscan(&vecsuf) 

dpsv[sf] [t] [c]=vecsuf 
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case 's' /* sigma a, a */ 

nsi— > nsi+1 printf(nsi) 

fscanf(Mrl,&udl, &spsufl) 

fscanf(&lr2,&ud2, &spsuf2) 

f lrl=0,lr2=l (case'901') sin -> sin+1, printf 

[ lrl=l,lr2=0 (case'910') sibarn — > sibarn+1, printf 

si[sf][t][c][lrl][udl]=spsufl 

si[sf][t][c][lr2][ud2]=spsuf2 

fscanf(&vecsuf) 

siv[sf ] [t] [c] =vecsuf 



case 'B' /* der Scalar dA, dA */ 

nB^ nB+1, dif— >dif+l printf (nB) 
fscanf(&lr, fevecsuf) 

lr=0 (dA) Bn -> Bn+1, printf(Bn) 
lr=l (<9A*) Bbarn — > Bban+1, printf (Bbarn) 
B[sf][t][c][lr]=vecsuf, printf 



case 'A' /* Scalar A, A */ 
nA— > nA+1, printf (n A) 
fscanf(&lr, &existno) 

lr=0 (A) An — >■ An+1, printf (An) 
lr=l (A*) Abarn — > Abarn+1, printf (Abarn) 
A[sf][t][c][lr]=existno, printf 



case 'C /* der de Scalar ddA, ddA */ 

nC-> nC+1, dif— > dif+2, printf(nC) 
fscanf(&lr, &existno) 

lr=0 (ddA) Cn -> Cn+1, printf(Cn) 
lr=l (<9cL4*) Cbarn — > Cbarn+1, printf(Cbarn) 
C[sf] [t] [c] [lr]=existno, printf 
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} /* c-running> */ 
} /* t-running> */ 
} /* sf-running> */ 

/* Expanding the product of superfields */ 
T[0]=0,1,- • • ,noterm[0J-l: Term No of Superfield 
T[1]=0,1,- • • ,noterm[l]-l: Term No of Superfield 1 

T[sfieldno-l]=0,l,- • • ,noterm[sfieldno-l]-l: Term No of Superfield sfieldno-1 



For every set (T[0],T[1],- • • ,T [sfieldno-1]), do the following box. 
printf (T[0],T[1],- ■ • ,T [sfieldno-1]) 

termscombineQ; 

thn=THnuml[0]; thbarn=THnuml[l]; 

if(thn >= 3) { printf("NO of theta: MORE THAN 2"); goto lab41; } 
if(thbarn >= 3) { printf("NO of theta-bar: MORE THAN 2"); goto lab41; } 

SORTOUTthBth(); 
thBthBthth(); 
SigmaContraction() ; 
VecSufChange(); 
lab41:; 



/* thn: No of 9 within a term; thbarn: No of 6 within a term; */ 
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10.2 termscombine 

printf("*** TERMSCOMBINE ***"); 
[InitializationJ 
c=0,l,- • • ,maxcompno2-l 
typel[c]=V 
thl[c][][]=99 
sil[c][][]=99 
dpsvl[c]=99 
psil[c][][]=99 
dpsl[c][][]=99 
Bl[c][ ]=99 
Al[c][ ]=99 
Fl[c][ ]=99 
Cl[c][ ]=99 



(T(O)term of sf=0 SuperField)x 
(T(l)term of sf=l SuperFieldjx 



(T(sfieldno-l)term of sf=sfieldno-l SuperField) 

— > One Term Combined 

Indices of a combined term 
"Gathered Term Check" 

No of Components, weight 1[0], weightl[l], DIFnuml, THnuml[ ], SIGnum[ 
], PSInuml[ ], Bnum[ ], Anuml[ ], Fnuml[ ], Cnuml[ 



10.3 SORTOUTthBth 

Moving th(9) , thbar(#), si(cr) and sibar(a) to the biggining part of a term. In 
the process, we obtain the sign change due to the Grassman algebra. 

PlusMinus : Sign change 

Here we use the function: idno(TP). 

idno(TP); a function which produces FieldID number correponding to a 
type TP. 
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10.4 thBthBthth 

First we change the data form of th[ ] [ ]. From : 

th l c JliJljJ= a 7 ij=0 or 1 
c=0,l,- • ■ ,nth-l 

to an another form: 

I[cJ=i, J [cj =j , Suf[cJ=a 

In this part, the program classifies terms by the power of 9 and 9. 

(1) indipendent of 9 and 9 (nth=0) 

(2) 9 a ,9 a ,9 6ll 9" (nth=l) 

(3) 9 a 9^9 a 9^--- (nth=2) 

(4) 9 a 9^ 9 a 9p, ■■■ ■ 9«9~P, 0*0^, ■ ■ ■ ; 9 2 ; 9 2 (nth=2) 

(5) 9 2 (9^9*),9 2 (9 a ,9 a ) (nth=3) 

(6) 9 a 9 f3 9^, ■ • • ; A 0^0 7 , ■ ■ ■ (nth=3) 

(7) 9 a 9^9^9 5 , • • • ; 9 & 9^9^9 S , (nth=4) 

(8) 9 a 9 $ 9~<9 A/ , • • • ; 9«9 /3 9 A <9 5 , (nth=4) 

(9) 9 a 9 /3 9^9*, • • ■ ; 9«¥9^9 & , (nth=4) 

The formula presented in SecJH] is exploited here. As the output, we ob- 
tain 

[OuTPut] 

weight, PlusMinus, Sign, Nthth, NthBthB, Nhalf 
This OutPut means the following factor 

(weight) * (_i)PlusMinus+Sign ^ ^ Nhalf ^ ^ Nthth ^ NthBthB 

appears in the process of this subsection, 'weight' is the starting one: 
weight=weight[0]+i weight [1]. Here we use SufChange(). 

SufChange(): Chiral suffix contraction using 5@, e Q/3 , e Q/ g. It is done for com- 
ponent c with c>nth or c>nth+nsi. 



10.5 SigmaContraction 

This routine is composed of 5 parts, (A),(B),(C),(D) and (E). 

(A) Change of the data form 
First we change the data form: (N=nsi-1) 
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a m °- 




rT m N a N 

a N 


si[nthjl2J[2j 




si[nth+lJ[2Jl2J 




,1 , '-inn 

silnth+nsi-lj[2j[2j 


siv[nth] 




siv[nth+l] 




r i i , • -i 1 

siv[nth+nsi-lj 


to an another form: 






c=0 


c=l 


c=N 


CsufflOj[OJ=a 




Csuff[lj[OJ=ai 




Csuff[NJlOj=a 7V 


Csuff[0][l]=d 




Csuff[l][l]=di 




Csuff[N][l]=djv 


Vsuff[0]=m 




Vsuff[l]=mi 




Vsuff[N]= m7V 


UD[0][0]=l(down) 




UD[l][0]=0(up) 




UD[N][0]=l(down) 


UDfOl [11=1 (down) 




UDfll[ll=Ofup) 




UD[Nlfll=Ofup) 



(B) Chiral-suffix-pair search 
Search for same suffixes in {Csuff[c][0]|c=0,l,- • • ,N}(odd number suffixes) 
and in {Csuff[c][l]|c=0,l,- • • ,N}(even number suffixes). When Csuff[cl][0]=Csuff[c2][0], 
then we call (cl,c2) chiral-suffix-pair (or, simply, Cpair). For the odd number 
case, we list them as 



CpairO[cpairnoO 


=oJLoj= 


cl 


CpairO[cpairnoO 


=o][i]= 


c2 (cl<c2) 


CpairO[cpairnoO 


=0][2]= 


=«o (common chiral suffix) 


CpairO[cpairnoO 


=1][0] = 


=c3 


CpairO[cpairnoO 


=!][!] = 


=c4 (c3<c4) 


CpairO[cpairnoO 


=1][2] = 




CpairO [NcpairO- 


l][0]=c 


* 


CpairO [NcpairO- 


l][l]=c 


+ (c*<c+) 


CpairO [NcpairO- 


a NcpairO~l 



For the even number case, 



CpairE[cpairnoE 


=oJLoj 


=clO 


CpairE[cpairnoE 


=o][i] 


=cll (clO<cll) 


CpairE[cpairnoE 


=0][2] 


=do (common anti-chiral suffix) 


CpairE[cpairnoE 


=1][0] 


=cl2 


CpairE[cpairnoE 


=1][1] 


=cl3 (cl2<cl3) 


CpairE[cpairnoE 


=1][2] 


=di 


CpairE[NcpairE- 


1][0] = 


c** 


CpairE[NcpairE- 


!][!]= 


C++ (c**<c++) 


CpairE[NcpairE- 


1][2] = 


O^NcpairE— 1 
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(C) Grouping a's 

A group of a's is defined to be a set of cr's which are connected by chiral- 
suffix or anti-chiral-suffix contractions. The number of groups within a term 
is assigned to be GrNum. This an important Graph index. The number of 
a's which each group has is stored in SigmaN[GR] (GR=0,1,- • • ,GrNum-l). 
Each group consits of some cr's specified their component numbers. They are 
stored in Group [GR][nseries] (nseries=0,l,-,SigmaN[GR]-l). In this process 
of grouping, the program traces the CpairO and CpairE defined in (B). Each 
group is determined by a series of cpairnoO's and cpainoE's. They are stored 
in CPseries[GR][ ]. 

(D) Vector-suffix-pair search 

Search the same vector-suffix in {Vsuff[c]|c=0,l,- • • ,N}. 
Vpair[0]lOJ=cl, Vpair[0]llJ=c2 
Vpair[l][0]=c3, Vpair[l][l]=c4 

Vpair[vpairno-I][0]=c*, Vpair[vpairno-l][l]=c+ 



(E) Classifying cr's part by indices (nsi, vpairno, NcpairO, NcpairE, closed 
chiral-loop) 

Using graph indices defined in the Sec|H we can classify cr's part. See the 
TABLE 2,3 and 4 in the text. Each class (we call here the finally classified 
place "class") has the number SigGraN which is the set of the graph indices. 
In each class, we reduce the expression using the graphical relations explained 
in SecEl Here the following important data are obtained and transfered to 
VecSufChange(). 

BranchN: Generally each class decomposes to some branches after the 
reduction using the graphical relations. We specify the number of branches. 
NoChangeN[br]: For each branch(br), NoChangeN is assigned, for the 
case that cr's part does not change, 1 for the case that cr's part changes. 
RedSigN[br]: For each branch, the decrease-number of a's after the reduc- 
tion is stored. This number is important for adjusting the output form, 
cstart: Normally this is fixed as cstart=nth+nsi 

LorenContN[br]: For each branch, specify the number of Lorenz contrac- 
tion, in other words, the number oft]™ 111 . 

SUFF[br][10][2]: Specify the vector suffixes appearing in 77's. 
MultiFac[br][2]: For each branch, a multiplicative factor is specified. 
silBR[br][c][2][2], sivlBR[br][c]: For each branch, specify the a's part. 
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epl[c][4]: For each branch, specify the totally-antisymmetric tensor when 
this term appears. 



10.6 VecSufChange 

Vector suffix contraction is peformed using rj mn rf s ■ ■ ■ . It is done for compo- 
nent c (c>cstart=nth+nsi). 

10.6.1 FinalOutPut 

11 Appendix B: Input and Output Examples 
11.1 Wess-Zumino Model 

We demonstrate the calculation of where $ is the chiral superfield, 

and obtain the component form. The next input data can be read from the 
graphical eqauations © and (J2J). 

INPUT DATA 

2 
6 






-1 


4 


t 


1 


s 


1 1 


112 51 


t 10 2 B 1 51 




1 





5 


t 


1 


t 


1 1 


t 112 


t 10 2 C 11 




2 





2 


t 


1 1 2 


P 


1 2 











1 


5 


t 


1 1 2 


t 


1 2 


t 1 


s 11 114 51 


q 1 


1 





3 


t 


1 1 2 


t 


1 2 


F 1 1 






1 





1 


A 


1 1 












6 























1 


4 


t 


5 


s 


1 5 


116 52 


t 10 6 B 52 




1 





5 


t 


5 


t 


1 5 


t 116 


t 10 6 C 1 




2 





2 


t 


5 


p 


1 5 











-1 


5 


t 


5 


t 


1 5 


q 7 


52 s 17 116 


52 t 


1 





3 


t 


5 


t 


1 5 


F 1 






1 





1 


A 


1 













OUT PUT 
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T[0]=0 T[1]=0 



****** TERMSCOMBINE **** 
th*th*thbar*thbar-term 

****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= l+i(0) PlusMinus= Sign=2 Nthth=l NthBthB=l Nhalf=2 

****** SigmaContraction **** 
SigGraN=201199 FinalOutPut: MultiFac=-2 + i(0) type2[c=6]= B B2[c=6,l]=51 
type2[c=7]= B B2[c=7,0]=51 

T[0]=1 T[l]=5 

****** TERMSCOMBINE **** 
th*th*thbar*thbar-term 

****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= l+i(0) PlusMinus= Sign=0 Nthth=l NthBthB=l Nhalf=0 

****** SigmaContraction **** 
SigGraN=0 FinalOutPut: MultiFac=l + i(0) type2[c=4]= C C2[c=4,l]=l 
type2[c=5]= A A2[c=5,0]=l 

T[0]=2 T[l]=3 

****** TERMSCOMBINE **** 
th*th*thbar*thbar-term 

****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= 0+i(-2) PlusMinus= 4 Sign=0 Nthth=l NthBthB=l Nhalf=l 

****** SigmaContraction **** 

SigGraN=l FinalOutPut: MultiFac=l + i(0) type2[c=4]= s si2[c=4,0,l]=7 

si2[c=4,l,l]=2 siv2[c=4]=52 type2[c=5]= p psi2[c=5,l,0]=2 type2[c=6]= 
q dps2[c=6,0,0]=7 dpsv2[c=6]=52 
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T[0]=3 T[l]=2 



****** TERMSCOMBINE **** 
th*th*thbar*thbar-term 

****** SORTOUTthBth **** lab50 at thBthBthth, Final result 
weight= 0+i(2) PlusMinus= 1 Sign=l Nthth=l NthBthB=l Nhalf=l 

****** SigmaContr action **** 

SigGraN=l FinalOutPut: MultiFac=l + i(0) type2[c=4]= s si2[c=4,0,l]=l 

si2[c=4,l,l]=4 siv2[c=4]=51 type2[c=5]= q dps2[c=5,l,0]=4 dpsv2[c=5]=51 
type2[c=6]= p psi2[c=6,0,0]=l 

T[0]=4 T[l]=4 

****** TERMSCOMBINE **** 
th*th*thbar*thbar-term 

****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= l+i(0) PlusMinus= Sign=0 Nthth=l NthBthB=l Nhalf=0 

****** SigmaContraction **** 
SigGraN=0 FinalOutPut: MultiFac=l + i(0) type2[c=4]= F F2[c=4,l]=l 
type2[c=5]= F F2[c=5,0]=l 

T[0]=5 T[l]=l 

****** TERMSCOMBINE **** 
th*th*thbar*thbar-term 

****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= l+i(0) PlusMinus=0 Sign=0 Nthth=l NthBthB=l Nhalf=0 

****** SigmaContraction **** 
SigGraN=0 FinalOutPut: MultiFac=l + i(0) type2[c=4]= A A2[c=4,l]=l 
type2[c=5]= C C2[c=5,0]=l 
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Gathering all terms, we obtain 

$t$| e2 - 2 == _l dmA * d ™ A (o 0) + }_Q2 A * . A (1? 5) + X _ A *tf A ( 5) !) 

-ix(-l) n(/°X/ (2,3)+<x(-l) V^/^XX (3,2) + F*F(4,4)(10) 

This is the Wess-Zumino Lagrangian. We donot ignore the total divergence 
here. 

11.2 Super QED 

The calculation of W a W a gives the kinetic terms of the photon and the 
photino in Super QED. The following input data can be read from the graph- 
ical eqauations (JHJ) and (JSJ). 

INPUT DATA 

2 
4 

0-1 1 1 Oil 

1 2 t 1 1 D 1 

0-1 4 s Oil 112 51 s 003 102 52 t 013 w 52 

51 

10 4 t 003 t 013 s 011 112 51 m 102 51 
4 

0-1 1 1 1 

1 2 t 1 D 1 

0-1 4 s 001 114 53 s 005 104 54 t 015 w 54 

53 

10 4 t 005 t 015 s 001 114 53 m 104 53 
OUT PUT 

T[0]=0 T[l]=3 

****** TERMSCOMBINE **** 
th*th-term 

****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 
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weight= 0+i(-l) PlusMinus= 2 Sign=0 Nthth=l NthBthB=0 Nhalf=0 



****** SigmaContraction **** 

SigGraN=l FinalOutPut: MultiFac=l + i(0) type2[c=2]= s si2[c=2,0,0]=l 

si2[c=2,l,l]=4 siv2[c=2]=53 type2[c=3]= 1 la2[c=3,0,l]=l type2[c=4] = 
m dl2[c=4,l,0]=4 dlv2[c=4]=53 



T[0]=1 T[l]=l 

****** TERMSCOMBINE **** 
th*th-term 



****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= l+i(0) PlusMinus= Sign=l Nthth=l NthBthB=0 Nhalf=0 



>£. >£. >£. 



SigmaContraction **** 
SigGraN=0 FinalOutPut: MultiFac=l + i(0) type2[c=2]= D D2[c=2]=l 
type2[c=3]= D D2[c=3]=l 



T[0]=1 T[l]=2 

****** TERMSCOMBINE **** 
th*th-term 



****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= 0+i(-l) PlusMinus=0 Sign=0 Nthth=l NthBthB=0 Nhalf=l 

****** SigmaContraction **** 
SigGraN=201199 FinalOutPut: MultiFac=2 + i(0) type2[c=4]= D D2[c=4]=l 
type2[c=5]= w dv2[c=5]=53 dvv2[c=5]=53 

T[0]=2 T[l]=l 

****** TERMSCOMBINE **** 
th*th-term 



****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= 0+i(-l) PlusMinus= Sign=l Nthth=l NthBthB=0 Nhalf=l 
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****** SigmaContraction **** 
SigGraN=201199 FinalOutPut: MultiFac=-2 + i(0) type2[c=4]=w dv2[c=4]=51 
dvv2[c=4]=51 type2[c=5]= D D2[c=5]=l 



T[0]=2 T[l]=2 

****** TERMSCOMBINE **** 
th*th-term 



****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= -l+i(0) PlusMinus=0 Sign=0 Nthth=l NthBthB=0 Nhalf=l 

****** SigmaContraction **** 
SigGraN=402299 
Branch No 

FinalOutPut: MultiFac=-2 + i(0) type2[c=6]= w dv2[c=6]=51 dvv2[c=6]=51 
type2[c=7] = w dv2[c=7]=54 dvv2[c=7]=54 
Branch No 1 

FinalOutPut: MultiFac=2 + i(0) type2[c=6]= w dv2[c=6]=52 dvv2[c=6]=51 
type2[c=7]= w dv2[c=7]=51 dvv2[c=7]=52 
Branch No 2 

FinalOutPut: MultiFac=-2 + i(0) type2[c=6]= w dv2[c=6]=52 dvv2[c=6]=51 
type2[c=7]= w dv2[c=7]=52 dvv2[c=7]=51 
Branch No 3 

FinalOutPut: MultiFac=0 + i(2) type2[c=5]= e ep(51, 52,54,53) type2[c=6] = 
w dv2[c=6]=52 dvv2[c=6]=51 type2[c=7]= w dv2[c=7]=54 dvv2[c=7]=53 

T[0]=3 T[1]=0 

TERMSCOMBINE **** 



^ ^ ^ ^ ^ ^ 



th*th-term 



****** SORTOUTthBth **** 
lab50 at thBthBthth, Final result 

weight= 0+i(-l) PlusMinus= Sign=0 Nthth=l NthBthB=0 Nhalf=0 

****** SigmaContraction **** 
SigGraN=l FinalOutPut: MultiFac=l + i(0) type2[c=2]= s si2[c=2,0,l]=l 
si2[c=2,l,l]=2 siv2[c=2]=51 type2[c=3]= m dl2[c=3,l,0]=2 dlv2[c=3]=51 
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type2[c=4]= 1 la2[c=4,0,0]=l 



Gathering all terms, we obtain 



/ N 



W a W a = 2i \y G XX {(0,3), (3,0)} + (1,2) + (2,1) 

- 1 - ^0 (2, 2)br0 - X -v mn 2 (2, 2)brl - l -v m 2 (2, 2)br2 - 2ie m > n > s > l ~v mn ■ \ sl (2, 2)br3 

-D 2 (1,1).(H) 
Adding WW, we obtain the Lagrangian as 

1 



C = -{-W a W a \ e * + W & W c 



=-\vmj+ i -( x/ a xy - x/^xx )+^ 2 - (12) 

where we do not ignore the total divergence. This is the kinetic term of the 
photon and the photino. 
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